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The effects of X-ray photoionization and heating on the 
structure of circumstellar discs 
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ABSTRACT 

We present the results of a theoretical study investigating the effects of photoioniza- 
tion and heating by X-rays on discs around low-mass stars. In particular we address 
the question of whether or not X-rays can drive a disc wind. First, we construct a 
1-dimensional "quasi-hydrostatic" model, which solves for the vertical structure in- 
troduced by X-ray heating. We consider uniform X-ray illumination of the disc, but 
the X-ray fluxes required to heat the disc significantly are much greater than those 
seen by recent observations. When the model is extended to consider heating from a 
central X-ray source we find that the 1-dimensional model is only valid very close to 
the star. We extend our analysis to consider a simple 2-dimensional model, treating 
the disc as a two-layered structure and solving for its density profile self-consistently. 
For T Tauri stars we are able to set a crude upper limit on the mass-loss rate that 
can be driven by X-ray photoevaporation, with a value of ~ 10 _13 g cm~ 2 s _1 . Our 
model is designed to maximise this value, and most likely over-estimates it signifi- 
cantly. However we still find a mass-loss rate which is less than that found in studies 
of ultraviolet photoevaporation. We conclude that in the presence of a significant UV 
field, X-ray driven disc winds are unlikely to play a significant role in the evolution of 
discs around low-mass stars. 
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1 INTRODUCTION 

The existence of discs around young low-mass stars is 
now well established and remains an important area of 
study, as disc evolution provides strong constraints on the- 
ories of both star and planet formation. The majority 
of young stars a re observed to have circumstellar discs 
at ages ~ 10 6 yr !|Sj i ri2rn_et alJ LL989]; iKenvon fc Hartmannl 
Il99.4 lHaisch. Lada fc Ladal |20Q1). These discs are opti- 
cally thick in the optical and infrared, an d typically have 
mass e s of a few percent of a s olar mass ijBeckwith et alJ 
Il99fl: lEisner fc Carpenter! 120031 ) . However these discs are 
not observed in most systems at ages of ~ 10 7 yr: gen- 
erally ordy_low-rnass_Mebris discs" remai n by this point 
(e.g. iMannines fc Sargent! Il997l : IWvatt. Dent, fc Greaves! 
2003). Thus disc lifetimes are expected to be ~ 10 6 yr. 

However few objects are observed to have prop- 
erties intermediate between those of disc-bearing, clas- 
sical T Tauri stars and discless, weak-lined T Tauri 
stars llKenvon fc Hartmannlll995l : lArmitage. Clarke fc Tout! 
1999). Thus very few objects are "caught in the act" 
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of dispersing their discs, and th erefore the dispersal 
timescale must be short, ~ 10 5 y r ( Simo n fc Pratd 1199,4 
IWolk fc Waited Il996t iDuvert et all l2000l) . The common 
mechanisms invoked as means of disc dispersal, such as 
magnetospheric clearing or viscous evolution, fail to sat - 
isfy this two-time-scale constraint JArmitage et al J 119991. 
Inst ead they produce pow er-law declines in disc proper- 
ties (Hartmann et al. 1998), thus always producing dispersal 
time-scales of the same order as the disc lifetime. 

One promising avenue of investigation, however, has 
been using photoevaporative heating as a means of driving 
disc dispersal. Models of the flows driven from steady discs 
by ultraviolet (UV) p hotoevaporation were first constructed 
by iHo llcnbach et alJ i 19941 , and have since been improved 
upon by several subsequent studies l^ grke^^Wclz 1996; 
iRicbling fc Yo7ke1IT997l: iKessel. Yorke fc Richliimlllflflrt: 



199c 



also the review by Hol lenbach. Yorke fc Johnstone! 120001. 
These models are characterised by mass loss outside some 
critical radius (usually referred to as the "gravitational ra- 
dius"), beyond which ionized material is heated to above 
the escape te mperature. The mass lo ss rate declines sharply 
with radius iJHollcnbach et al.lll994l find an R~ 5 ' 2 scaling 
beyond the gravitational radius), and so most of the mass 
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loss o ccurs at or very n ear to the gravitational radius. Re- 
cently |Fbnt_etjil] J2004F) have performed detailed hydrody- 
namic modelling of this process, and whilst they find the 
mass-loss profile around the gravitational radius to be mod- 
ified somewhat, their results are in broad agreement with 
those from the earlier models. 

Initially these models failed to reproduce observed 
disc lifetimes, as the low mass-loss rates meant that 
the disc was dispersed much too slowly. However 
more re cent studies dClark e. Gend rin fc Sotomavorl 
|22^]j JMatsuvama. Johnstone fc Hartmannl I2H 



lArmitaee. Clark e & Palla 2003) have coupled the pho- 
toevaporative mass loss with viscous evolution, and found 
that a number of observed disc properties can be satisfied 
by these mo d els. In particular, the "UV switch" model of 
IClarke et alJ (|2001) was able to reproduce the two-time- 
scale behaviour described, and so these models provide an 
attractive area of study. However a number of problems still 
exist. The models fail to satisfy observational constraints on 
the sub-mm emissio n associated with t he outer parts of the 
disc at late times l|Clarke et alJ 120011. and detailed SED 
modelling has not yet been performed. Further, and more 
importantly, the details of the Lyman continuum emission 
from the central star, emission required to drive the pho- 
toevaporation, are poorly understood. E ven the magnitude 
of the emission is poorly constrained JGahm et al.1 Il979t 



llmhoff fc Appenzellerlll987ll . In a previous paper we have 
shown that emission from the accretion shock is unlikely to 
drive disc photoevaporation (Alexander. Clark e fc Pringld 
2004), and it is not known if stellar chromospheres can 
provide the requisite flux of ionizing photons. Further, given 
the problems associated with observing at wavelengths 
immediately shortward of the Lyman limit, it is not clear 
when, or if, this issue will be resolved. 

A logical step beyond UV photoionization models is to 
look at the effect of X-ray ionization and heating on disc 
evolution. X-rays are obviously capable of ionizing disc ma- 
terial and, unlike in the UV case, the X-ray emission from 
young stellar objects and T Tauri stars (her eafter TTS) 
has been well studied (e.g. reviews by iNeuhauserl Il997t 
iFeigelson fc Montmerlelll999l and references within) . X-rays 
are emitted from both magnetic reconnection events near 
to the star and also from protostellar outflows/jets, in both 
cases with luminosities potentially large enough to influence 
disc evolution. Further, recent observations of young stel- 
lar clusters containing massive stars (of type earlier than 
05) have found diffuse X- ray emission throughout the clus- 
ters IITownslev et al . 2003). The effect of this external X-ray 
flux on disc structure and evolution may also warrant inves- 
tigation. 

Previous studies of X-ray/disc interactions 
n 



(e.g. lOlassgold. Naiita fc Igeal 119971: lleea fc Glassgoldl 
Il999l: iFromang. Terauem fc Balbusl |2002) have focused 
mainly on using X-rays as a means of sustaining the low 
levels of io nization required to drive the magnetorotational 
instability (Ba lbus fc Hawlevlll991l) . or studied the effect 
of heating to low temperatures (< 1000K) on the observed 
spectrum (Glasseold & Naiita 2 0011 : IGorti fc Hollenbachl 
1200411 . Here we make a preliminary study of X-ray heating 
as a means of driving a disc wind. The structure of the paper 
is as follows. In Section [5] we review relevant observations 
and X-ray physics, and then set up our basic disc model. In 



Section |3 we use this 1-dimensional (1-D) model to study a 
disc subject to uniform X-ray illumination, and then apply 
the model to illumination from a central source in Section 
[I] In this case the model is of limited validity, and so in 
Section |S| we extend it to a simple 2-dimensional (2-D) 
model. In Section HJ we discuss the caveats and limitations 
that apply to the model, and in Section |7| we summarize 
our conclusions. 



2 DISC WINDS AND X-RAY HEATING 

2.1 Background 

Previous studies of photoevaporatively driven disc winds 
have focused on the case of UV photoionization, where 
the disc is heated by photons with energies of 13.6eV. 
This has bee n studied in detail b y a number of 
authors (e.g. iHollenbach et all Il994l : iRichling fc Yorkd 
W97|_jFon^^alJ|2004j; see also the review chapter by 
Hollenbach et al. 2000). The 2-D radiative transfer scheme 
of Hollenbach et al. (1994) finds that a central UV source 



produces a hot ionized atmosphere above the surface of the 
isothermal disc. They find that the diffuse field (ionizing 
photons produced by recombinations of Hi atoms in the at- 
mosphere) dominates over the direct field from the central 
source. Beyond a gravitational radius, the thermal energy of 
the particles in the atmosphere becomes greater than their 
gravitational potential energy, and so beyond this gravita- 
tional radius ionized particles escape from the disc: this is 
photoevaporative mass-loss. 

By comparison to the UV case the case of X-ray heating 
is more complex, due to the more complicated nature of the 
interaction of X-rays with matter. The physical processes 
involved in X-ray absorption by diff use gas are described 
in detail bvlKro lik & Kallmar j (I1983J) (see also the review s 
bv lFeigelson fc Montmerlelll999l and blassgold et alJl200d1 . 
Here we merely summarise the salient points. 

The X-ray emission around Young Stellar Objects 
(YSOs) and TTS is well studied, and is due to three 
distinct mechanisms. First, and best-studied, is the X- 
ray emission from the central star and its immediate sur- 
roundings. This is thought to arise from magnetic field 
reconnection events in the stellar corona and magneto- 
sphere, and typically produces an X-ray spectrum consis- 
tent with thermal bremsstra hlung (e.g. iNeuhauseil Il997l : 
IFeigelson fc Montmerld [.1999) . This emission peaks at at 
around 2keV, and YSOs and TTS typically exhibit steady 
X-ray luminosities of around 10 28 -10 30 erg s . They also 
exhibit rapid variability and occasional "flaring", when 
large changes in X-ray luminosity, sometimes several or- 
ders of magni tude, occur on timesca les of order hours 
(e.g. llmanishi. Kovama fc Tsu boi 2001). 

The second source of X-ray emission around young stars 
is from the jets and shocks observed in protostellar outflows. 
Such sources are visible at other wavelengths as Herbig-Haro 
(HH) objects, and produce X-ray emission with similar prop- 
erties (s pectrum, luminosity etc.) to that produced by the 
corona iFeigelson fc Montmerld Il999l) . However these jets 
and shocks are far from the central star and disc, at dis- 
tances typically tens to hundreds of AU. This geometric di- 
lution of the energy, combined with probable extinction due 
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to absorbing material between the HH shock and the disc, 
means that this X-ray emission is unlikely to have any sig- 
nificant effect on the disc. Therefore we do not include this 
emission in our modelling. 

The third mech anism is a very recent discovery. 
iTownslev et alJ l|2003r) observed X-ray emission throughout 
M17 and the Rosette Nebula, and attributed it to thermal 
bremsstrahlung emission from the diffuse ionized ISM. They 
observed integrated luminosities of order 10 33 erg s _1 from 
areas on the sky of order ~ 50pc 2 . Little variation was ob- 
served across the spatial extent of the clusters, and the ob- 
served spectrum was slightly softer than that seen from mag- 
netospheric X-rays (peaking at 0.8-1.5keV). However this 
emission has only been in observed in clusters containing 
very massive stars. Indeed t he lack of any similar emisssion 
in the Orion Nebula Cluster JFeigelson et all 20031 indicates 
that stars of spectral type 05 or earlier may be required in 
order to produce this emission in young stellar clusters. It 
follows that this cannot be a universal mechanism for disc 
dispersal, as low-mass star forming regions such as Taurus- 
Auriga contain no massive stars. It may be relevant in some 
cases, however, and is also much simpler to analyse than the 
case of emission from a central source. 

When soft X-ray emission is incident on the gas in a 
disc incident X-ray photons are usually absorbed by heavy 
elements, such as O, C or Fe. A single photoelectron is pro- 
duced, which then collides with neutral material in the gas, 
generating typically 30 "secondary" electrons. Heavy ele- 
ments can also undergo the Auger effect, in which the ex- 
cited ion produced by the initial photoabsorption undergoes 
2-electron decay, and direct (collisional) charge exchange be- 
tween neutral hydrogen or helium and ionized heavy species 
can also occur. A key difference bewteen the X-ray and UV 
cases is the details of the incident radiation field. Whilst in 
the case of UV ionization the recombination (diffuse) field 
dominates, the probability of an X-ray photon being emit- 
ted by a recombination is small JGlasseold et alll997n . Thus 
only the direct field is important in X-ray heating, and the 
heating is local to the initial X-ray absorption. Consequently 
care must be taken with regard to the geometry of the prob- 
lem. Of the 3 cases discussed above, we neglect the emission 
from the HH outflow. First, we consider the geometrically 
simpler case of unifor m illumination from a d iffuse plasma, 
of the type observed bv lTownslev et ah! ( 2003) . Having stud- 
ied the effects of uniform illumination we then discuss the 
more promising, but also more complicated, case of X-rays 
generated by magnetic activity near to the surface of the 
central object. 

Throughout the following calculations we adopt an in- 
cident X-ray spectrum consistent with that of 10 7 K opti- 
cally thin bremsstrahlung, with a peak at ~ 0.7keV. These 
energies are typical of both the diffuse and coronal emis- 
sion observed (see discussion above) . In order to ensure that 
the models consider only the effect of the X-rays and not 
the "tail" of the spectrum at lower energies, we subject the 
bremsstrahlung spectrum to an exponential cutoff at ener- 
gies E < O.lkeV. 



2.2 Basic disc model 

As a first iteration we make use of a simple hydrostatic disc 
model. We set up a steady disc, and then consider the ef- 



fects on the structure of X-ray heating. Initially we treat the 
radiation as being vertically incident on the disc, and so we 
make the approximation that the disc can be treated as a 
series of concentric, non-interacting annuli. Here we define 
our steady disc model. 

The surface density, E(-R), of the steady disc is given 
by 

M d 



E(R) 



2tt R s R 



exp(—R/R s 



(1) 



This represents a power-law decline in surface density sur- 
face density (c.f. lBeckwith et aljll990l iBell et al.lll997fl pro- 
portional to 1/-R, tapered exponentially at around a scaling 
radius R a . We note however that our results do not depend 
strongly on the exact form of E(7?). We also note that this is 
similar to the form of theoretically derived time-dependent 
solut i ons for the surface density (e.g-lLynden-Bell fe Pringla 
11974 iHartmann et al.lll99Sl: IClarke et alJl200lfl . We adopt 
a typical disc mass of Md = 0.01M Q . We adopt a scaling 
radius of R s = 10AU initially, and later consider the effect 
of varying this parameter. 

At a given radi us R the verti cal structure is that of an 
isothermal disc (see |Pringlelll98lI) . The general equation of 
hydrostatic equilibrium perpendicular to the plane of the 
disc is 

i op _ d f gm. \ (] 

pdz dz\ ^Ri + z 2 ) { ' 

However for small vertical displacements z above or below 
the midplane we can approximate the gravitational potential 
to be harmonic, and so the equation of hydrostatic equilib- 
rium simplifies to 

ldP _ GM,z 
~p~dz~ ~ W 

As shown in IPringld II1981T) we can solve this to find the 
isothermal vertical density structure, which is a Gaussian 
distribution in z: 

■2 



(for z < R) 



(3) 



p(z) = p (R)exp - 



2H 2 (R) 
Here the density po{R) is given by 

fm gCg) 

po{R) = T^m 

where the scale height H(R) is defined to be 

?3\ 1 > 2 



H(R) = 



cj{R)R- 
GM* 



and c s is the local isothermal sound speed, 

P _ kT(R) 
p pniH 



c 2 s (R) 



(4) 



(5) 



(6) 



(7) 



Thus choosing the disc temperature to be a function of R 
determines the disc structure uniquely. The behaviour of 
the disc temperature with radius is not completely under- 
stood, and so power-laws are usually adopted. Simple repro- 
cessing of stellar radiation by a thin disc results in a T ex 
R~ 3 ' 4 scaling \ Adams fc Shu 1986), and accretion-powered 
"self-luminosity" results in the same scaling relationship 
JLvnden-Bell fe Pringlenl974) . However such a relationship 
does not fit observed data, and tend s to under-predict the 
flux observed at long wavelengths IIKenvon fe Hartmannl 
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119871) . Much work has been done on this (e.g. see discus - 
sions in lKenvon fc Hartm ann 1987t iHartmann et alJ ll998). 
and the consensus is that a T oc R~ ' scaling law is both 
consistent with observational data and theoretically justi- 
fia ble (for example the "flared reprocessing disc" proposed 
bvlKenvon fc HartmanrJll987l) . Therefore, again following 
IClarke et alj (120011) . we adopt a power-law scaling of the 
form 



T(R) = 



-1/2 



T 



(8) 



with the normalisatio n condition To = 100K at Ro = 1AU 
iJBeckwith et alJll99Ci) . We adopt a fiducial stellar mass of 
M„ = IMq throughout. 



3 UNIFORM ILLUMINATION 

3.1 Modelling X-ray heating 

In order to study what we expect to happen to such a disc 
when it is subjected to, and subsequently heated by, ionizing 
radiation, we have first mad e use of the CLOUDY photoion- 
ization code JFerlandlll996l) . Given a fixed, static density 
structure this code solves the equations of thermal and ion- 
ization equilibium to find a unique solution. Only a density 
structure, incident flux spectrum and chemical compsition 
must be specified as input parameters 1 . When run, CLOUDY 
returns temperature and ionization profiles for the heated 
region. It also highlights the main heating and cooling pro- 
cesses. 

Initially, we treat the disc as a series of concentric, non- 
interacting annuli heated vertically from above. This reduces 
the problem to a series of independent 1-D problems. We 
take the vertical density profile of the steady disc and use 
this as an input to the CLOUDY code. If the sound crossing 
timescale, given by t s ~ c s /H{R), is much shorter than the 
thermal equilibrium timescale then we can assume that the 
disc will "relax" to a non-isothermal vertical structure much 
faster than the heating processes change. Thus we can take 
the temperature profile obtained from CLOUDY and solve the 
equation of hydrostatic equilibrium in the vertical direction 
to find a new, updated density profile. This is in turn used 
as an input to CLOUDY, and so we iterate towards a self- 
consistent solution. In practice the sound crossing timescale 
is always found to be at least 10 times less than the thermal 
equilibrium timescale, and so this approximation is valid. 

When we solve the equation of hydrostatic equilibrium 
(Equation |3J in the non-isothermal case, we find that the 
particle number density n(z) is given by 



, •, Tr 

n(z) =n R —— 



T(z 



■ exp 



m H GM, r fj,(z')z' 



i? s 



1 ~R 



Tt 



■dz 



(9) 



for z > zr. In this equation T(z) is the temperature, n(z) the 
particle number density, fi(z) the mean molecular weight, 
and the subscripted values Tr and ur refer to the values of 
these functions evaluated at some reference point zr, which 
is introduced as a boundary condition. This solution is pre- 
sented in full in Appendix 1X1 

In practice, as can be seen from Equation [5] the choice 

1 We assume solar abundances throughout. 



of initial density structure has a strong effect on this pro- 
cess, and not all initial structures converge towards a self- 
consistent solution. In order to do this we adopt a "two- 
component" initial profile, consisting of a cold central core 
(that of the unperturbed, isothermal disc) and a hot "at- 
mosphere" in the region where we expect the disc to be 
heated. The core has the standard isothermal structure 
above (Equation 0J. The heated atmosphere has a similar 
structure, but this time with a scale height 

/ C 2 fl 3 \ 1/2 _/ fcT hot fl 3 \ 1/2 
Hhot [gM-J U"»hGmJ (10) 

For the regions of interest, temperatures of 1000-10, 000K 
are typical, and so we adopt a value of T^ ot = 3000K for 
this initial guess. We note however that as long as the pro- 
cedure converges the results do not depend on this value. 
For the inital profile we define a reference height zt at 
which temperature changes from the midplane temperature 
to Thot = 3000K. We maintain pressure equilibrium at this 
point, and so the inital density profile shows a contact dis- 
continuity at this point. The criterion for determining zt is 
discussed below. 

Our model is limited by the fact that the CLOUDY code is 
only valid for temperatures > 3000K, and so we cannot solve 
explicitly for the structure all the way down to the disc mid- 
plane. Consequently we are forced to extrapolate from the 
height at which the temperature falls to 3000K down to the 
point at which it reaches the midplane temperature (zr in 
Equation^. We extrapolate linearly, keeping the gradient of 
T(z) constant from the stopping point of the code (3000K) 
down to the midplane. In doing this we neglect the exis- 
tence of an extended column of "warm" material (~ 1000K). 
Other studies, while exploring somewhat different parameter 
spaces, have found that such an extended column probably 
does result when circumstellar discs are heated by X-rays 
jGlasseold fc Naiital200ltlGorti fc Hollenbachl2004l) . which 
may modify the resulting mass- loss rates somewhat. 

In this initial model we evaluate the location of the 
transition point empirically, fixing zt so that the T(z) pro- 
file can be extrapolated smoothly. (In Section |K| we define 
an explicit criterion for the location of the transition point, 
but this is not necessary at this stage.) We always define 
zr (Equation^ to be the point at which the density profile 
first diverges from that of the isothermal midplane. An ex- 
ample of the iterative process is shown in Fig0 Typically 
this process take 5-10 iterations to converge. 



3.2 Results 

We have first considered uniform illumination of the type 
expected in the massive star clusters described in Section 
12.11 We evaluated the disc model described above with a 
uniform X-ray flux Fx — 10~ 2 erg s _1 cm ~ 2 . The vertical 
structure of the disc was evaluated at a number of different 
radii spanning all of the expected "static" region, as well as 
two models which are outside the gravitational radius. We 
evaluated the model at radii of IORq (= 7.0 x 10 11 cm), 2OR.0 
(= 1.4 x 10 12 cm), 5.0 x 10 12 cm, 1.0 x 10 13 cm, 5.0 x 10 13 cm, 
1.0 x 10 14 cm, 2.0 x 10 14 cm, 3.0 x 10 14 cm, 4.0 x 10 14 cm. At 
each radius we truncated the disc at the point where the 
density of the initial profile fell to 1cm -3 . 

The density and temperature profiles evaluated by the 
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terations Example: Central Illumination at R=1QR 




z / cm 

Figure 1. Example of iterations towards final density profile. 
The inital profile is shown as a dashed line, the subsequent (7) 
iterations as dotted lines, and the final profile as a solid line. 



F x =10 erg s 




Figure 2. Density (heavy line) and temperature (light line) struc- 
tures evaluated at R = 5 X 10 12 cm. The solid line represents the 
region calculated explicitly by the model, the dashed line the re- 
gion where we extrapolate linearly in T(z), and the dotted line 
the isothermal midplane. 



model at R = 5 x 10 12 cm are shown in Fig|3 (These are typ- 
ical of the results at all values of R.) The temperature rises 
smoothly from the midplane temperature to around 6,000- 
10,000K. The gas cannot be heated to temperatures much 
greater than this without being collisionally ionized, and this 
only happens in regions where the particle number density is 
very low (< 10 cm -3 ). Thus the region of interest is mostly 
neutral, with ionization fractions of a few percent at most. 
By far the dominant heating process in this region is photo- 
electronic heating, which accounts for > 90% of the observed 
heating. Near to the disc surface direct charge exchange is 
also significant, typically accounting for a few percent of the 
total heating. The contribution of Auger electrons to the 
total heating is negligible, and the X-rays produced by re- 
combinations (mostly Fe K«) are negligible compared to the 
incident flux. At the densities of interest here the dominant 
cooling mechanism is metal line emission. Only at densities 



> 10 10 cm -3 do 2-body cooling processes become significant, 
so we do not see significant 2-body cooling here. 

Due to the large column density of the region heated by 
X-rays, some three orders of magnitude greater than that for 
UV photons, the heated region covers a large spatial extent, 
and in fact expands to a height where the z <C R approxi- 
mation (Equation |3J is no longer valid. This is discussed in 
Section |S| below. Also, we note that the outer three models 
are evaluated outside the scaling radius of 10AU (Equation 
0, in the region where the surface density is falling expo- 
nentially. However, as long as E(-R) is large enough that the 
unperturbed disc is extremely optically thick at a given ra- 
dius R, the structure of the cold core has little effect on the 
structure of the heated region. This condition is satisfied 
throughout the disc, and so the choice of scaling radius does 
not affect our results significantly. 

Whilst this model is only valid in the inner, static re- 
gion, a similar structure would be expected in the outer, 
flow region (although it would obviously be altered some- 
what by the flow of heated gas away from the disc). Thus we 
can discuss the potential form of mass-loss in a qualitative 
manner, even if our model cannot account for the process 
in detail. As the heating occurs over a much larger column 
than in the UV case, with a gradual decrease in temperature 
with decreasing height above the midplane, we would expect 
the mass-loss profile to differ from the UV case. In the UV 
case, where the disc is divided sharply into hot and cold re- 
gions at the ionization front, most of the mass is lost at, or 
very close to, the gravitational radius. However the smoother 
temperature profile seen in the our model means that com- 
parable mass loss rates exist at a number of different radii, 
due to the fact that the escape temperature decreases as 
1/R. As a result, the point at which the thermal energy of 
the gas becomes greater than its gravitational energy occurs 
at a different height above the midplane at each different 
radius, and so the mass-loss need not be concentrated at 
the gravitational radius. Indeed, evaluating pc s at the point 
where T(z) = T osc for the outermost three models (which 
lie at the 8000K gravitational radius and beyond) results in 
very similar values. We find mass-loss rates per unit area 2 
of 8.6 x 10 _17 g cm -2 s" 1 at R = 2.0 x 10 14 cm, 2.0 x 10~ 16 g 
cm -2 s _1 at -R = 3.0 x 10 14 cm and 2.5 x 10~ 16 g cm -2 s" 1 
at R = 4.0 x 10 14 cm. However, as mentioned above, this 
is really over-interpreting our explicitly hydrostatic model, 
and so these numbers are not considered to be accurate. 

In the very low density regions of the disc, at large z, 
most of the gas is ionized. Where the gas is mostly ion- 
ized, large X-ray fluxes can heat the gas to ~ 10 6 K, poten- 
tially providing some mass-loss at radii as small as 0.1 AU. 
However the densities in these regions are so low that the 
maximum mass-loss rates (per unit area) attainable are 2-3 
orders of magnitude lower than those at at 10AU, so this 
effect can safely be neglected. 

In principle this seems to be a promising avenue of in- 
vestigation, with a real potential to influence disc evolu- 
tion. However, the problem with these models is the lack 
of X-ray photons in real systems. Our models show that 



2 

a mass 



For reference, the UV models of Hollcnba ch et al.l i 1994) derive 
dss rate per unit area of 3.88 x 10 — 13 g cm -2 s — 1 at 
cm. 
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significant heating on ly occurs for va l ues of Fx > 10 4 erg 
s _1 cm ~ 2 . However. JTownslev et al.l 1120031) observe fluxes 
of ~ 2 x 10~ 7 erg s" 1 cm ~ 2 in M17 and ~ 4 x 10" 8 erg 
s _1 cm ~ 2 in the Rosette nebula. These values are some 
three orders of magnitude too small to affect the disc struc- 
ture significantly, and so it seems unlikely that such diffuse 
bremsstrahlung emission will play a significant role in disc 
evolution. 



4 ILLUMINATION BY A CENTRAL SOURCE: 
1-D MODEL 

We now consider the case of a disc illuminated by X-rays 
from a central source. The X-rays emitted from YSOs and 
TTS have been observed for many years, and a wealth of 
observational data exist. As discussed above, the X-rays are 
subject to rapid variability on small scales, and also occa- 
sional "flaring" on scales of several orders of magnitude. If 
we neglect flaring, however, it is a reasonable approximation 
to treat the source as having a quasi-steady X-ray luminosity 
of 10 28 -10 30 erg s _1 . Here we assume that the central source 
has a luminosity of Lx = 10 30 erg s _1 and radiates isotrop- 
ically. Whilst this may not be valid, the detailed geometry 
of the X-ray emission is not well enough understood to jus- 
tify any other selection. The X-rays are thought to originate 
from the reconnection of field lin es in the upper parts o f the 
magnetosphere, and so, following iGlasseold et al.l al997lt . we 
place the X-ray source at R = 0, z s — IORq. The validity 
of this appproximation is discussed in Section |S| below. 

In the case of optical heating of discs by a central 
source fe.g.lAda ms fc Shull986|:|jKenyon fc Hartm annll987t 
IChiang fc Goldreichlll997l : lDullemond et alj|200ll) it is com- 
monly assumed that photons are absorbed at the cylindrical 
radius where they hit the disc. Disc heating then occurs 
through a diffuse radiation field whose net direction, given 
a disc-like geometry, is approximately vertically downwards. 
Thus the flux at each radius is given by: 



L x =10""ergs 



F x = 



4tt7? 2 V# 2 + if 



(11) 



where the first factor in this equation represents simple geo- 
metric dilution of the energy and the second term is the ver- 
tical component of the incident flux, given a source height 
z a — IORq. The procedure then is to compute the ver- 
tical hydrostatic equilibrium structure at each cylindrical 
radius, subject to this boundary condition on the incident 
flux. Although we follow this procedure here as a first guess 
at the structure of an X-ray irradiated disc, we note that 
it is less valid in the X-ray case. The reason for this is 
that th e diffuse radiation fie ld is negligible in the case of 
X-rays (IGlassgold et alJ ll997). so that heating instead pro- 
ceeds purely by attenuation (due to absorption) of the in- 
cident X-ray beam. Obviously a self-consistent treatment 
would involve an iterative solution for the 2-D structure of 
the disc, subject to hydrostatic equilibrium along each ver- 
tical slice and thermal equilibrium along each ray path from 
the source. However as a first step we adopt this 1-D ap- 
proximation. 

The expression for Fx in Equation 1111 assumes that the 
X-ray flux is downwardly incident on the disc surface. Thus 
it is valid as long as the disc scale height is much less than 




Figure 3. As Fig|2] but for the central illumination model- 
evaluated at R = IORq. The change in gradient of T(z) at 
z ~ 6 X 10 10 cm is due to the increased rate of 2-body cooling 
as the density increases. 

Spatial variation of disc density 



I 






R / R G 

Figure 4. Plot of the density distribution produced by the 1-D 
model. The greyscale strips are the density profiles evaluated by 
our model, with dotted contours added as a guide to the eye. 
Contours are drawn at n = 10 16 , 10 15 , 10 14 . . . 10 8 cm- 3 . The 
position of the X-ray source is marked by circles. 



z s . If the disc expands to a height greater than z 3 this ap- 
proximation breaks down. Thus we can use the 1-D model 
described in Section^lto solve for the disc structure close to 
the star, subject to this condition on the validity of Equation 



EH 



4.1 Results 

We evaluate these models near to the central source, at radii 
of IORq, 2OR , 30R Q , 40Rq and 5OR . In general, the 
structures are very similar to those seen in the uniform ra- 
diation field modelled in Section |3] above. The heating rate 
is much greater, as the local X-ray fluxes so close to the 
source are much higher (~ 10 4 erg s _1 cm ~ 2 ), but the heat- 
ing and cooling processes are very similar. Due to the high 
X-ray fluxes we were forced to truncate the initial density 
profiles at the point where the density falls to 10 4 cm~ 3 to 
avoid numerical problems in solving for the ionization bal- 
ance. The temperature and density profiles resulting from 
the R — IORq model are shown in Fig[3] The steepening 
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Two— temperature fit at R=1OR 



z / cm 

Figure 5. Two-temperature fit to the R = lOR© model. The 
solid line is the density profile evaluated by the model. The dashed 
line is the best-fitting two-temperature profile, with parameters 
z T = 4.21 X 10 10 cm and T hot = 6200K. The contact disconti- 
nuity in the two-temperature density profile at zj- arises because 
pressure equilibrium is maintained over the temperature discon- 
tinuity. 



of the T(z) curve as z falls to ~ 6 x 10 10 cm is due to the 
increased significance of 2-body cooling as the density rises. 
The overall effect on the structure of the inner disc is shown 
in Fig[l 



4.2 Two-component fit 

As seen in Fig^] at radii > 3OR.0 the disc expands to a 
height comparable to z s and so the approximation of down- 
ward incidence is no longer valid. Also, we note that the 
column density along the line-of-sight from the source to zr 
at R = 5OR.0 is ~ 10 26 cm~ 2 . This corresponds to an optical 
depth of ~ 10 4 at IkeV, and so it is clear that attenuation 
of the X-ray flux through the disc will be significant at large 
radii. Thus in order to evaluate the global structure of a disc 
heated by a central X-ray source we must both consider the 
problem in two dimensions and also account for attenuation 
of the flux through the disc. 

As a first step along the path to a simple 2-D model we 
note that the vertical structures at these radii can be rel- 
atively well-fit by a "two-component" profile, parametrized 
in the same way as the initial profile described in Section 
13.11 (with free parameters zt and Th ot ). Through a simple 
least-squares fit we found that a best-fitting temperature of 
6200K provided a good fit to the structures at both 10R© 
and 2OR.0 (see Fig|SJ. We will make use of this fact in the 
next section, where we consider a simple 2-D model. 



5 ILLUMINATION BY A CENTRAL SOURCE: 
SIMPLE 2-D MODEL 

In order to study the behaviour of the disc at larger radii, 
in particular nearer to the gravitational radius, we seek to 
construct a "toy" 2-D model of the disc structure. As seen 
in Sections [3] and 2] above, the vertical density structure 



of the disc is reasonably well approximated by a two-layer 
structure, with the cold disc beneath a layer heated by the 
X-rays. We have already seen that the vertical structure can 
be reasonably well- fit by a two-temperature structure, and 
so an approximate structure of this form provides a good 
approximation in the parts of the disc which control the 
photoevaporative mass-loss rate. The key issue controlling 
the mass- loss rate is the location of the transition point zt, 
which marks the bottom of the heated region: to first order, 
the mass- loss rate per unit area varies as pc s , where p is the 
density at zt- Thus we now seek to solve for the density 
n(R, z) in a self-consistent manner. 

The location of the transition point zt is determined 
by the depth to which the X-ray penetrate the disc. In the 
optically thin limit, with little or no absorption of the inci- 
dent flux, the importance of X-ray effects is determined by 
the parameter 



nd? 



(12) 



where Lx is the X-ray luminosity, d the distance from 
the source and n the particle number density at dis- 
tance d. In essence this parameter is a measure of 
the photon-to-particle density ratio at any given point 
(e.g. iTarter. Tucker fc Salpeterl Il969t iGlasseold fc Naiital 
2001). Thus in our model X-ray heating is significant where 
£ is greater than some critical value, and where £ is less than 
this value it is not. 

Complications arise, however, when optical depth ef- 
fects become important, as the term Lx/d 2 in Equation 1121 
is effectively a photon flux and so is reduced by extinction 
effects. At a single photon energy the optical depth is sim- 
ply given by N/N c , where N is the column density along 
the line-of-sight from the source and N c is the critical col- 
umn density for which the optical depth is unity. However 
as the integrated X-ray photoionization cross-section varies 
strongly with photon energy E, so does N c . For solar abun- 
dances, N c varies as 



N c = 



4.4 x 10 21 cm" 2 



(13) 



— V 

ikevy 

with a be st-fitting value of a = 2.485 IjGlassgold et alJ 
ll997U200dV Thus, with a varying incident continuum (10 7 K 
bremsstrahlung) and a varying absorption cross-section 
we must integrate over photon energy to find the total 
attenu ation. Using the formulation of iKrolik fc Kallmanl 
( 1983 ), updated to us e the improved cross-section fits of 
lOlassgold et, al.l l|1997l) . the attenuation factor A is given 
by the integral: 



/"OO 

Jh(r, Xo) — I x~ a exp (— x — rx~ a ) 

J in 



dx 



(14) 



where xo = Eo/kTx is the lower cutoff of the spectrum and 
r here is the optical depth to photons with energies at the 
peak of the continuum (0.7keV). The first term represents 
the falling absorption cross-section, and the second is a com- 
bination of the exponential fact ors in the the bremsstra hlung 
continuum and the extinction. lOlassgold et alJ l|1997l) show 
that Jh can be well-approximated by the expression 



Jh(r) = At a exp ( —Bt 



with best fitting parameters of A 



0.800, 



(15) 



0.570, 
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X— ray heated disc 



Disc structure 
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Figure 6. Number density at the base of the heated region, , 
as evaluated by the 2-D model, plotted as a function of radius. 
The solid line shows a simple R~ 2 powerlaw. The only significant 
deviation from this line is at small radii. At larger radii the atten- 
uation factor J/j is approximately constant with radius, implying 
that most of the attenuation occurs close to the star. The slight 
deviation from the R~ 2 powerlaw at very large radii is due to the 
effects of the scaling radius R s on the density profile, and is not 
significant. 



B = 1.821 and 6 = 0.287 (for solar abundances and keV 
energies). Consequently, in the case where attenuation ef- 
fects are significant the location of zt is not given by a 
critical value of £, but rather a critical value of the modified 
parameter 



e = 



LxJh 

nd 2 



(16) 



Thus our procedure for evaluating the disc struture 
is as follows. Firstly we interpolate the inner structure, 
evaluated in Section [I] above, onto a regular (R, z) grid 
over the range R = 1O-2OR0. We adopt grid spacings of 
2R = 1.39 x 10 11 cm in R and 1 x 10 9 cm in z. In this region 
the approximation in Equation 1111 is valid, and so we use 
the results of the 1-D model to evaluate the critical value of 
£' at the transition point to be £' c — 4.8 x 10 -8 erg cm s . 
This uncertainty in this value is approximately 10%. The 
attenuation through the disc to a radius Ro depends only 
on the structure of the disc at R < Ro- Consequently, if we 
know the structure in the region R < Ro we can evaluate Jh 
by integrating the column along the line-of-sight to (Ro, z), 
and using the fit in Equation ll5l In this manner we evaluate 
all values of £'(z) at the radius Ro. The height above the 
midplane which gives £'(z) = £.' c is adopted as the location 
of zt at radius Ro, and we then create the resulting vertical 
structure at Ro using the two-component profile described 
above (with the best-fitting temperature of Th ot = 6200K). 
We can then move to the next grid cell and repeat the pro- 
cess, and thus iterate outward though the disc to find the 
global disc structure. 



5.1 Results 

We evaluate the model described above out to a radius 
of 50AU. The procedure finds values of £' c (zt) accurate to 




R / R Q 

Figure 7. Structure of the X-ray heated disc evaluated by the 
2-D model. Density contours are drawn at n = 10 15 , 10 14 , 
10 13 . . . 10 3 cm -3 . The location of the X-ray source at R = 0, 
z s = IORq is marked by circles. 



within 5% of the critical value in all cases, and within 1% 
for all radii greater than 10 13 cm= 0.67AU. This translates 
to an uncertainty of less than 10% in the value of the den- 
sity at this point. Runs at higher spatial resolution (x2, x5) 
over a smaller radial extent do not change the numerical ac- 
curacy significantly, and so we consider the procedure to be 
numerically stable. 

The value of n(R,zr) (i.e. the number density at the 
base of the heated region) as a function of radius R is shown 
if Fig|S] and the global disc structure is shown in FigQ We 
see from Figlol that n(R, zt) diverges from a R~ 2 power- 
law at radii < IOORq. Looking at Equation 1161 we see that 
this implies that the attenuation factor A is constant at 
radii > IOORq. This in turn implies that almost all of the 
attenuation of the X-ray flux occurs very near to the star. 
The slight divergence from the R~ 2 power-law at large radii 
is a numerical effect due to the value of the scale radius R s 
and is not significant. 

As before, a first-order estimate of the mass-loss rate 
per unit area can be made by evaluating pc s at the tran- 
sition point. At 2.85 x 10 14 cm= 19.1AU, the gravitational 
radius for 6,200K material, this gives a value of 2.6 x 10 _13 g 
cm~ 2 s _1 . Larger disc scale radii R s can reduce this some- 
what, but it remains constant to within a factor of 3 for 
R s < 50 AU, and larger scale radii still make little or no 
difference. Due to the manner in which we determine zt, 
neglecting the possibilty of an extended column of "warm" 
material (see Section f3, 11 . this value of the mass-loss rate is 
regarded as an upper limit. As the disc density falls off as 
exp(—z 2 /2H 2 ) even small changes in the position of zt can 
result in large changes in the value of n at this point, with 
correspondingly large changes in the mass-loss rate. Thus we 
regard this value as an upper limit, with the true mass- loss 
rate being sm aller, possibly significantly so. The mass-loss 
rate derived by Hollcnbach ct al. ( 1994) for UV photoevapo- 
ration is 3.88 x 10 _13 g cm -2 s _1 at the gravitational radius. 
Thus we find that the mass-loss rate of an X-ray driven disc 
wind is as best comparable to the UV rate, and more likely 
considerably less significant. 
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As a consistency check we took the density structure 
along different lines-of-sight from the source and used these 
as inputs for the CLOUDY code. This test indicates that the 
location of the transition point is accurate to within 10% at 
radii less than 1AU. At larger radii the test confirms that 
we place the transition point at too low a vertical height. 



6 DISCUSSION 

There are obviously a number of simplifications in these 
models. In the 1-D model there is firstly the treatment of 
the problem as a series of concentric, non-interacting annuli. 
This is obviously unphysical and in reality we would expect 
some radial heat transfer, almost certainly outward through 
the disc. However in the limit of vertical heating this radial 
heat transfer will be negligible by comparison, and so this 
approximation is sound. 

In addition, we point out that our modelling is explic- 
itly hydrostatic, yet attempts to deal with the hydrodynamic 
behaviour of the disc. The treatment of the static region, 
inside the gravitational radius, is valid, as the sound cross- 
ing timescale is always much less than the thermal equilib- 
rium timescale. However outside the gravitational radius the 
model is clearly not valid, as material will flow away from 
the disc at a speed comparable to the local sound speed. 
Thus the calculations at radii > 2.0 x 10 14 cm are suspect, 
and are considered to be illustrative only. 

In the case of a central X-ray source, its location is not 
well-established. The X-ray emission is thought to arise from 
magnetic reconnection events in the stellar magnetosphere, 
and so should be located in the upper region of the magne- 
tosphere. Previous studies hav e used either a point source 
located at -R = 0, z — 10R© JGlasseold et al.lll 997t). or a 
ring at R = 5-lORp , z = 5-10R^ ftgea fc Glassgoldlll999l; 
iFromang et ali .2002'). We adopt the point source for simplic- 
ity, and do not investigate the effect of moving the source. 
However it is unlikely to affect the results significantly. We 
find that the heated disc causes significant attenuation of 
the X-ray flux for all R > 3OR0. Moving the source closer 
to the disc will likely increase the heating at small radii, but 
the effect far from the source will be small. 

There are several other, more minor problems with the 
1-D model which also merit discussion. One further com- 
plication is that the X-ray-heated discs expand significantly 
in the vertical direction, often well-past the point to which 
the z <C R approximation (Equation [3J ceases to be valid. 
In essence our disc resides in a 1-D harmonic potential, and 
thus we overestimate the potential at large z. This will lead 
to the model somewhat underestimating the vertical expan- 
sion of the disc. However the densities of the disc in this 
region are very small, and so little mass is affected by this 
simplification. Also, the manner in which the vertical expan- 
sion is evaluated does not explicitly conserve mass. However 
the increase in mass due to this simplification only occurs at 
a level corresponding to around 10 -8 of the total disc mass, 
and so we neglect its effect. 

There is also the issue of grains. Most of the X-rays are 
absorbed by heavy elements, which in our model are dis- 
tibuted evenly throughout the disc (i.e. with constant grain- 
to-gas ratio). However a significant fraction of the heavy 
elements can be "locked up" in grains, and it has been pro- 



posed that g rains sediment towards the disc midpl ane (see 
the review bv lBeckwith. Henning fc Nakagawal2000l and ref- 
erences within). This has the knock-on effect of moving the 
X-ray absorption closer to the disc midplane, thus possibly 
changing the distribution of heating in the vertical direction. 
This issue is not addressed by our model. 

Our 2-D model is very simplistic, and is only really a 
"toy" model of the process. Our criterion for the transi- 
tion point is correct, in as much as it correctly identifies the 
point at which the X-rays can heat the disc to temperatures 
> 3000K. However, as discussed in Section 13.11 above, the 
way in which we use the location of this point to create a 
vertical density structure is somewhat idealised. By assum- 
ing that the temperature continues to fall smoothly we im- 
plicitly assume that there is no extended column of "warm" 
material at ~ 1000K. If such a column exists then the point 
at which the temperature is high enough to drive mass-loss 
is moved to higher z, and therefore somewhat lower density. 
Also, the nature of the two-temperature fit is such that we 
tend to under-estimate the density at z>zt- The approx- 
imation of a single temperature in the heated region means 
that we under-estimate the density where the temperature 
is much greater than our best-fitting Th ot = 6200K. Con- 
sequently we under-estimate the extinction along the lin- 
of-sight to the outer parts of the disc, and therefore over- 
estimate the depth to which the X-rays penetrate. A more 
realistic simulation would incorporate varying temperatures 
along the different lines-of-sight, but this is beyond the scope 
of our simple model. We note however, that both of these 
simplifications have been made in such a way as to max- 
imise the mass-loss rate. In addition, we have neglected the 
possibility that there is any attenuation of the X-ray flux 
between the source and the disc (apart from that arising in 
the disc atmosphere), for example by an outflow or an accre- 
tion column. Any reduction in the photon flux reaching the 
disc would obviously reduce the mass-loss rate still further. 
We also note that we have adopted a high X-ray luminosity 
of Lx — 10 30 erg s _1 . As can be seen from Equation 1 161 the 
density at the base of the heated region, and therefore the 
mass-loss rate, essentially scales as n oc Lx- Thus a smaller 
X-ray luminosity, more typical of TTS, will reduce the mass- 
loss rate still further. Moreover this linear drop-off is much 
faster than the square-root dependence o n the ionizing flux 
found in models of UV photoevaporation (iHollenbach et all 
1 19941 . 

Consequently we consider our mass-loss rate to be a 
robust upper limit to the mass-loss that can be driven by 
X-ray heating, with the true value being lower, possibly con- 
siderably so. In the even of a significant Lyman continuum 
flux being present X-ray heating effects will be negligible. 
However in the absence of a strong Lyman continuum the 
X-rays may be able to provide mass-loss at a rate that is 
significant in disc evolution models. Further investigation 
of this problem is still needed, using models that can treat 
the thermal structure of X-ray irradiated gas over a large 
dynamic range of temperatures (i.e. 100-10, 000K). 

Comparisons to previous studies of X-ray/disc inter- 
actions are of limited use only, as t hese studies have had 
very different goals. Some studies fe.g.lGlassgold et aljl997l : 
llgea fc Glassgoldlll999l : IFromang et al]l2002ft have focused 
on sustaining the very low levels of ionization required to 
drive the magnetorotational instability. These studies find 
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that this low level of ionization is provided mostly by hard 
X-rays (3-10keV), which penetrate much further into the 
disc than the softer X-rays in which we are interested. Very 
few of these hard X-rays are absorbed near to the disc sur- 
face, but there are too few hard photons to provide signifi- 
cant heating. Further, the heated region resulting from our 
model is in fact optically thin to these hard photons, so 
the hard X-rays will be absorbed in t he cold region close 
to the midplane. More recent studies (Glassgol d fc Naiital 
l200ll:lGorti fc Hollenbacbll2004h have looked at the observa- 
tional signatures of X-ray heating. These arise at lower z and 
higher density than the region in which we are interested, 
and so again we are not able to compare to these directly. 
Thus our results are consistent with these previous studies, 
bearing in mind the important caveat that we are interested 
in very different physical effects. 



7 CONCLUSIONS 

We have modelled the effects of X-ray heating on the struc- 
ture of circumstellar discs, looking in particular at the ques- 
tion of whether X-rays can drive a disc wind. We first con- 
structed a 1-D model and used it to study the case of uni- 
form X-ray illumination. Such X-ray emission does exist in 
some regions of massive star formation, but at a level that 
is some 4-5 orders of magnitude too low to heat the disc 
significantly. We have then looked at heating from a cen- 
tral X-ray source, using parameters typical of T Tauri stars. 
Initially we used our 1-D model, but found that it breaks 
down very close to the central star. Thus we extended our 
analysis to consider a simple 2-D model in which X-rays 
reaching the outer disc are attenuated by passage through 
the X-ray heated atmosphere of the inner disc. We find that 
X-rays are significantly attentuated in this way, with most 
of the attenuation occuring within ~ lOOR© of the central 
star. The critical parameter in determining the mass-loss 
rate is the location of the transition point zt, which marks 
the lower boundary at which X-ray heating is significant. We 
have therefore constructed a simple 2-D model, solving for 
zt in a self-consistent manner. The model is constructed in 
such a way as to maximise the mass-loss rate, and most likely 
over-estimates the mass-loss significantly. However even our 
upper limit of ~ 10 _13 g cm -2 s _1 is less than that derived 
for ultraviolet photoevaporation. A more realistic 2-D sim- 
ulation, incorporating the thermal response of the gas over 
a wide range of temperatures, is necessary in order to inves- 
tigate this process fully. In the absence of a significant UV 
flux an X-ray driven disc wind may be significant. However 
if a significant UV flux is present, it seems likely that the 
UV-photoevaporation dominates. Given the constraints pro- 
vided by current observational data and the results of this 
work it seems unlikely that X-ray photoevaporation is a sig- 
nificant factor in the evolution of circumstellar discs around 
low-mass stars. 
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If we now integrate from some reference point zr to a height 
z > zr, we find 



log P + log 



m H 



m H GAL r ^'K dz > (A7) 



k R? J ZR T(z 



and so 



P(z) = P fl exp — - / -TjTT-r-dz (A8) 



k R 3 J ZR T(z') 

where the subscript r indicates values at z = zr. Finally 
we use the ideal gas law again to find the expression for the 
density profile as a function of the temperature profile and 
ionization structure. 

, , T R (ran GNU [ z p(z')z' , , N 

n(z) = riR — t N exp ' ' •'■ 



T(z) 



k R 3 



~M d >) (A9 » 



APPENDIX A: SOLVING THE EQUATION OF 
HYDROSTATIC EQUILIBRIUM 

The equation of hydrostatic equilibrium (Equation|3]m Sec- 
tion l2.2t governing the vertical structure of an accretion disc 
is 



GhUz 



ldP _ 

~p ~dz ~ R 3 

and the ideal gas law is: 



(for z < R) 



nkT 



pkT 
firriH 






m H 
k ' 



(Al) 



(A2) 



where fi is the mean molecular weight, mu is the mass of the 
hydrogen atom and k is Boltzmann's constant. Substutiting 
P from Eg uation I A2I gives 



dP 
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(pdT + Tdp 


pT dp 










dz 


m H 


\ p, dz p dz 


p? dz 



(A3) 



and so we can substitute this expression into Equation IA1I 
to find 



1 dT I dp 1 dp _ mu GM, pz 
T~dz~ + ~p~d~z~ ~Js,~dz~ ~ k RFY 



Rearranging, we see that 

d (. pT 
=> — log i— 
dz \ p 



m H GM„ pz 



k R 3 T 

and we can use the ideal gas law again to find that 

d ( run \ m H GM* pz 

^_(^ log P + log _j =_____ 
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